Strong field double ionization : The phase space perspective 
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We identify the phase-space structures that regulate atomic double ionization in strong ultrashort 
laser pulses. The emerging dynamical picture complements the recollision scenario by clarifying the 
distinct roles played by the recolliding and core electrons, and leads to verifiable predictions on the 
characteristic features of the "knee" , a hallmark of the nonsequential process. 

PACS numbers: 32.80.Rm, 05.45.Ac 



One of the most striking surprises of recent years in 
intense laser-matter interactions has come from multi- 
ple ionization by intense short laser pulses: Correlated 
(nonsequential) double ionization rates were found to be 
several orders of magnitude higher than the uncorrelated 
sequential mechanism allows. This discrepancy has made 
the characteristic "knee" shape in the double ionization 
yield versus intensity plot into one of the most dramatic 
manifestations of electron-electron correlation in nature. 
The precise mechanism that makes correlation so effec- 
tive is far from settled. Different scenarios have been 
proposed to explain the mechanism behind ionization [1, 

iiis s a s s m m m m p ^MM m 

and have been confronted with experiments [19j, [20( , the 
recollision scenario 0, H| , in which the ionized electron 
is hurled back at the ion core by the laser, being in best 
accord with experiments. In Fig. [TJ a typical double ion- 
ization probability as a function of the intensity of the 
laser field is plotted. Similar knees have been observed 
in experimental data Q, i, M, EL EL EL El US HI 
and successfully reproduced by quant al computations on 
atoms and molecules lEl. l27l l28l |29| . In a recent series of 
articles [1, HI, [H, EL 29T 3oT 31 1 characteristic features 
of double ionization were reproduced using classical tra- 
jectories and this success was ascribed to the paramount 
role of correlation |12|. Indeed, entirely classical inter- 
actions turn out to be adequate to generate the strong 
two-electron correlation needed for double ionization. 

In this Letter, we complement the well-known recol- 
lision scenario by identifying the organizing principles 
which explain the statistical properties of the classical 
trajectories such as ionization probabilities. In addition 
to the dynamical picture of the ionized electron provided 
by the recollision scenario, we connect the dynamics of 
the core electron and the energy flow leading to double 
ionization to relevant phase space structures (periodic or- 
bits or invariant tori). The resulting picture leads to two 
verifiable predictions for key points which make up the 
knee in Fig. [TJ namely the laser intensity where nonse- 
quential double ionization is maximal and the intensity 
where the double ionization is complete. 

We work with the classical Hamiltonian model of the 
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FIG. 1: Double ionization probability for Hamiltonian ([TJ) 
as a function of the intensity of the field I for u = 0.0584 
a.u. The vertical lines indicate (in green) the laser intensity 



4.57 x 10 W • cm where our dynamical analysis 



predicts the maximum of nonsequential double ionization, and 
(in blue) the intensity I c — 1.86 x 10 16 W ■ cm -2 where the 
double ionization is expected to be complete (see Fig. [4| . 



helium atom with soft Coulomb potentials [32|, [33| ■ The 
Hamiltonian is given by [101 ] : 



H(x,y,p x ,p y ,t) = f + ^- + ( 
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where x, y and p Xl p y are the positions and (canoni- 
cally conjugate) momenta of each electron respectively. 
The energy is initially fixed at the ground state £ g = 
—2.24 a.u. 34]. The laser field is modeled by a sinu- 
soidal pulse with an envelope, i.e. E(t) = Eq f(t) smut 
where Eq is the maximum amplitude and to the laser fre- 
quency chosen as u> — 0.0584 a.u. which corresponds to 
a wavelength of 780 nm. The pulse envelope fit) is cho- 
sen as a trapezoidal function with 2-4-2 laser pulse shape 
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(the ramp-up lasts two cycles, the plateau four, and the 
ramp-down two) [HI, [H, IS, US] • Typical ionizing trajec- 
tories of Hamiltonian fl} show two qualitatively different 
routes to double ionization : nonsequential double ion- 
ization (NSDI), where the two electrons leave the core 
(inner) region at about the same time, and sequential 
double ionization (SDI), where one electron leaves the 
inner region long time after the other one has ionized. 

We first analyze the dynamics of Hamiltonian |T]) with- 
out the field (Eq = 0) using linear stability properties 
such as obtained by the finite-time Lyapunov (FTL) ex- 
ponents [35[. FTL maps quantify the degree of chaos 
in phase space and highlight invariant structures. With 
each initial condition of a particular ensemble (here taken 
on the plane (x,p x ) with y = 0, the other coordinate 
p y being determined by the energy condition H = E g ), 
we associate the coefficient log |A(i)|/i where X(t) is the 
largest eigenvalue of the Jacobian obtained at time t from 
the tangent flow |35j. A typical FTL map is depicted in 
Fig. [2] for Hamiltonian ([I]) without the field. It clearly 
displays strong and global chaos by showing fine details 
of the stretching and folding of trajectories In par- 
ticular, there are no regular elliptic islands of stability 
contrary to what is common with Hamiltonian systems 
on a bounded energy manifold. By examining typical tra- 
jectories, we notice that the motion of the two electrons 
tracks, at different times, one of four hyperbolic periodic 
orbits, denoted O x .\, O x .2, O y> i and O yi 2 and displayed in 
Fig. [21 Their period is 29 a.u., i.e., much shorter than the 
duration of the laser pulse (of order 800 a.u.). These four 
orbits are symmetrical [reflecting the symmetries of the 
Hamiltonian |[T])] : The two periodic orbits X) i and O x ^2 
[outer projections in the plane (x,p x )] have on the plane 
(y,p y ), the same projections as the periodic orbits O y> \ 
and O y ^2 [inner projections in the plane (x,p x )]- Conse- 
quently, a typical two-electron trajectory is composed of 
one electron close to the nucleus (the "inner" electron) 
and another further away (the "outer" electron), with 
quick exchanges of the roles of each electron. We will see 
below that this distinction is crucial: Since the contribu- 
tion of the field-electron interaction to Hamiltonian (Q]) 
is proportional to the position, the action of the field is 
larger for the outer electron, while the inner electron is 
mainly driven by the interaction with the nucleus. 

Single ionization- By switching on the field, the outer 
electron is picked up and swept away from the nucleus. 
Consequently, its effective Hamiltonian is : 

2 

Hi = y + E Q xf(t) sinwt. (2) 

We notice that Hamiltonian Hi is integrable. Its solu- 
tions are approximately composed of linear escape from 
the nucleus (at time to) modulated by the action of the 
field [1, [H, [I?) (see the red trajectory in Fig. [3]). 

For the inner electron, the effective Hamiltonian con- 
tains the interaction with the nucleus and with the laser 
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FIG. 2: FTL maps of Hamiltonian {TJ without the field at 
time t = 43 a.u. in the plane (x,p x ) with y = 0. The projec- 
tion of some important periodic orbits (as continuous curves) 
and their respective Poincare sections (as dots) : O x ,i dashed 
dotted line (pink online), O x ,i dashed line (blue online), O y ,i 
dotted line (green online), O h ,2 full line (cyan online). 

field : 

H 2 = ^-- 2. + yE sin cot. (3) 

In the absence of the field (Eq = 0), TL2 is also inte- 
grable and the inner electron is confined on a periodic 
orbit. Since it stays close to the nucleus, its approximate 
period is 2tt/^/2 obtained from the harmonic approxima- 
tion, as observed in Fig. [3] 

Sequential double ionization (SDI)- Once an elec- 
tron has been ionized (usually during the ramp-up of 
the field), the other electron is left with the nucleus 
and the field. Its effective Hamiltonian is H.2- A con- 
tour plot of the electron excursions after two laser cycles 
and a Poincare section of H2 are depicted in Fig. Q] for 
J = 5 x 10 15 W-cm~ 2 . They clearly show two distinct re- 
gions : The first one is the core region which is composed 
of a collection of invariant tori which are slight deforma- 
tions of the ones obtained in the integrable case H.2 with- 
out the field. This elliptic region is organized around a 
main elliptic periodic orbit which has the same period as 
the field 2tt/u w 107.6 a.u. In this region, the electrons 
are fairly insensitive to the field, and do not ionize. The 
second region is the one outside the core where trajec- 
tories ionize quickly. It corresponds to sequential double 
ionization. In between these two regions, any weak inter- 
action (with the outer electron for instance) may move 
the inner electron confined on the outermost regular tori 
(but still inside the brown elliptical region) to the outer 
region where it ionizes quickly. 

If the laser intensity / is too small, then the phase 
space is filled with invariant tori and no sequential double 
ionization can occur because the motion is regular. The 
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FIG. 3: Two typical trajectories of Hamiltonian {T} for / = 
10 15 W ■ cm for initial conditions in the ground state energy 
of the helium atom. The two positions (x in red and y in 
blue) are plotted versus time (expressed in laser cycles). Note 
the vastly different vertical scales of the two panels. The 
recollision mechanism is seen in both panels : In the upper 
one, the recollision (at the end of the panel) brings in enough 
energy to ionize the inner electron. In the lower panel, the 
recollision energy insufficient to ionize the inner electron - the 
electrons exchange roles instead. 



sequential double ionization probability depends then on 
the size of the regular region around the elliptic periodic 
orbit, and hence on /. We have computed the location 
and the stability of this periodic orbit for ui — 0.0584 a.u. 
using a Newton-Raphson algorithm [35j . When it exists, 
this periodic orbit stays elliptic in the whole range of 
intensities we have considered. On the stroboscopic plot 
(with frequency lu) the periodic orbit is located at y = 0. 
In Fig. [3J the momentum p y of the periodic orbit on 
the stroboscopic plot is represented as a function of /. 
We notice that for a large set of intensities in the range 
[10 14 , 10 16 ]W-cm~ 2 , this periodic orbit is located close to 
p y = 0. For intensities larger than I c — 1.86 x 10 16 W • 
cm~ 2 , the periodic orbit does not exist, and no major 
islands of regularity remain. Therefore, it is expected 
that the sequential double ionization probability is equal 
to one in this range of intensities, as observed on the 
probability curve on Fig. [T] 

Nonsequential double ionization (NSDI)- As noted 
before, when the field is turned on, its action is concen- 
trated on only one electron, the outer one, as a first step. 
The field drives the outer electron away from the nucleus, 
leaving the inner electron nearly unaffected by the field 
because itsposition remains small. From the recollision 
process @, y| , the outer electron might come back close 
to the nucleus during the pulse plateau, if the field am- 
plitude is not too large. In this case, it transfers a part 
of its energy to the inner electron through the electron- 
electron interaction term. From then on, two outcomes 
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FIG. 4: Upper panel: Contour plot of the electron location 
y(t) at time t = 215.2 a.u. (2 laser cycles), of Hamilto- 
nian ((3]) for I = 5 X 10 15 W ■ cm -2 . Poincare sections (stro- 
boscopic plot) of selected trajectories in the elliptic central 
region are also depicted. The color code is on a logarithmic 
scale. The inset shows a projection of this periodic orbit at 
/ = 1.7 x 10 16 W • cm~ 2 in the (t/,p !/ )-plane. Lower panel : 
Momentum of the central periodic orbit (on the Poincare sec- 
tion) of Hamiltonian (|3j as a function of the laser intensity. 
The vertical line on the lower panel indicates the intensity 
Ic = 1.86 x 10 16 W ■ cm" 2 such that for / > I c , complete 
unhindered SDI is expected. 



are possible : If the energy brought in by the outer elec- 
tron is sufficient for the other electron to escape from the 
regular region (as in Fig. [31 upper panel), then it might 
ionize together with the outer electron. The maximum 
energy £ x of the outer electron when it returns to the in- 
ner region (after having left the inner region with a small 
momentum po close to zero) is obtained from Hamilto- 
nian ^ and is £ x — kU p , where U p = E /(4:Ui 2 ) is the 
ponderomotive energy and k = 3.17... is the maximum 
recollision kinetic energy in units of U p @, |H, H3] • We 
complement the recollision scenario (which focuses on the 
outer electron) by providing the phase space picture of 
the inner electron : In order to ionize the core electron, 
the energy brought back by the outer electron has to be 
of order of the energy difference between the core (y = 0) 
and the boundary of the stable region (y = y m ) of Jii 
(see Fig. SJ) which is equal to 
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A good approximation to y rn = y m (Eo) is given by the 
value where the potential is locally maximum, i.e. Eq — 
2y m /(yf n + l) 3 / 2 . The equal-sharing relation which links 
the classical picture of the outer electron x with the one 
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of the inner electron y, 

defines (through an implicit equation) the expected value 

(c) 

of the field Eq for maximal NSDI, because it describes 
the case when each outer electron brings back enough en- 
ergy to ionize the inner electron, while remaining ionized 
itself. However, fulfilling this energy requirement does 
not guarantee NSDI: The outcome depends on the num- 
ber and efficiency of recollisions. For uo — 0.0584, the 
predicted value of the amplitude E^f 1 as given by Eq. ([5]) 
corresponds to an intensity of 1^ = 4.57 x 10 14 W- cm -2 
which agrees very well with the simulations shown in 
Fig. [TJ In a wide range of frequencies, an accurate ex- 
pansion of Eq is obtained from Eqs. |l|)-((5]) and given 
by 



for sufficiently small to. To leading order the correspond- 
ing intensity varies as to 2 . For uj = 0.0584, the approxi- 
mate intensity given by Eq. © is 4.60 x 10 14 W • cm~ 2 
which is in excellent agreement with 1^°' . When the 
field Eq is too small, then the outer electron cannot gain 
enough energy to ionize the inner electron. When the 
field Eq is too large, then the outer electron does not rec- 
ollide since it leaves the interaction region nearly linearly. 
These two limits explain the bell shape of the resulting 
nonsequential double ionization probability, which, when 
put together with the monotonic rise of the SDI probabil- 
ity at higher intensities, adds up to the knee in question. 
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